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1. Introduction 

In the Tevatron collider and the LHC, multiple final states are observed frequently. In order to 
check the standard model and looking forward to finding new physics beyond the standard model 
|], ||, |3|, (|, we need to explore the problem of how to calculate the cross section efficiently and 
conveniently. In the past, the cross section is evaluated in the 3-dimensional momentum space § 
and people have developed quite mature numerical techniques. For the applications of programs 
Madgraph, Pythia, AlpGen and Sherpa, the reader can check references, for example, 0. 
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On the other hand, enormous progress have been made in the evaluation of one-loop amplitudes 
|| . One of such progress is the unitarity cut method originally proposed in || ||. With the 
twistor program initiated by Witten |1(J, the double cut phase space integration has been reduced 



to algebraic manipulation through the holomorphic anomily [|TT| , |12[ [U|, [14]]. Inspired by this 



simplification, in our first paper [|L5], we have explored how to apply the spinor integration method to 
the evaluation of the cross section for massless case. There are some obvious advantages comparing 
with the momentum integration method. First, the three-dimension momentum space integration 
can be reduced to just one-dimensional integration and furthermore for the massless case, the 
integration region is just [0, l] 2 . Secondly, in the calculation, every step is manifestly Lorentz 
invariant, thus we obtain compact analytic expressions. 

Continuing our study for the massless case, in this paper we focus on the massive case. We 
will see that if all the mass is set to zero, the massless case will be reproduced. Different from 
the massless case, the integration variable L is no longer a null momentum. So we can't apply the 
spinor method directly. However this problem have been solved in the unitarity cut method [|16|, [L7| . 
More accurately, we can write 

J d 4 L = J dz J d 4 £5 + (f)(2£- K); L = £ + zK, (1.1) 

where K is a fixed vector and z is a real number. Through this decomposition, we establish the 
general framework for massive case by the spinor method. 



In our first paper, we have emphasized the advantages of using the spinor method In the 
massive case, the constrained three-dimensional momentum space integration still can be reduced 
to an one-dimensional integration, plus possible Feynman integrations. In every step, we get a 
scalar type of integrations, which are Lorentz invariant. Furthermore, the integration region can be 
written directly. Though it's not simply [0, 1] like the massless case, it's only the simple functions 
of mass and energy. 

The outline of this paper is as follows. In section 2, we first briefly review the 4D unitarity cut 
method. Then we take the Faddeev-Popov trick to establish the general framework. 

In section 3, we apply our method to the pure phase space integration for two, three and four 
out-going particles as well as some simple examples to demonstrate the main idea and feature. 
These are the basis for practical and more complicated applications. 

In section 4, we calculate two practical examples and summarize some experience of performing 
the integrations. 

A summary of our results with some comments is given in section 5. 



2 The result of unitarity integration maybe written as one Feynman parameter integration over rational functions. 
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2. Framework to use spinor method 



In this section, we will setup the spinor integration method for massive particles in 4D. Then we 
apply this method to the phase space integration where the integration region (i.e, the j dx) of one 
dimensionless parameter is determined by the kinematical discussion. This region corresponds to 
the boundary of the whole phase space of outgoing momenta. One important difference, compared 
to the massless case, is that the integration region will be functions of masses of outgoing particles. 



2.1 The spinor integration method for massive cuts 

Here, we briefly review the spinor integration method for massive cuts (or sometimes called the 



"unitarity cut method") flq , |17| [18| . The Lorentz-invariant phase space (LIPS) of a massive double 
cut is defined by inserting two ^-functions representing the cut conditions: 



d 4 



(2.1) 



where I is the internal loop momentum and K the total momentum through the unitarity cut. 
Because I is a massive momentum, to use spinor integration method we need decompose it as 



+ zK, 



0; 



dH= / dzdH 5 + (f)(2£- K). 



(2.2) 



where i is a null 4-momentum, and can be expressed with spinor variables as 



tP xl , P x ~ x = AA; J dH S + (f ) = J (A A) A A J tdt. 



(2.3) 



Under this decomposition Eq. (|2.1|) becomes 



dzdH8 + (t)(2i ■ K)5 + (z 2 K 2 + 2zK ■ i - m{)5 + {{l - 2z)K 2 - 2K ■ i + mf - mf 2 ) 
dz ((1 - 2z)K 2 + m\- ml)5 + (z{\ - z)K 2 + z{m\ - m\) - m\) 

2 i ™2 _2 i^2 , ™2 ™2 



(A A) 



A A 



1 - 2z)K 2 + m{-mi _ (1 - 2z)K 2 + m{-mi. 



\\K\X 



MKW 



(2.4) 



In the first line of the result it depends only on the variable z, so we can use the ^-function to 
eliminate z as follows 



z± 



(K 2 + m\ — ml) ± a/A [K, m 1; m 2 ] 



2K 2 



(2.5) 
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where we have defined 



A[K, mi, 1712} = (K — m l — m 2 ) — 4m 1 m 2 . 



(2.6) 



Between the two solutions of z, only one should be taken. To see that, we make a kinematical 
analysis. Choose a center-of-mass frame such that 



K = {E> 0,0,0,0), 



(a, b, 0,0), K-£ = (E-a,-b,0,0). 



The mass-shell conditions require a 2 — b 2 = m 2 and (E — a) 2 — b 2 = m 2 ,, so a = (E 2 + m 2 — mfyjlE. 
In the decomposition i = I + zK, because the positive light cone with S + (£) have been chosen, we 
can write 

£=(|6|,6,0,0), I=(\b\ + zE,b,0,0) 



Then |6| + z.E = a. This means that only z_ is retained 3 . In the following of this paper, we always 
refer to z as 2_ , if it's not explicitly illustrated. 
Then Eq.(|2.4j) becomes 



(A A) 



A A 



1 - 2z)K 2 + m\-m 2 2 _ (1 - 2z)K 2 +m\-ml 



A LPT A 



A LPs: A 



(2.7) 



Eq. ( p.7|) is our final form for the spinor integration with massive double cuts. For convenience we 
define 



a[i^,mi;m 2 ] -/3[i^;mi,m 2 ] /3-PT 2 
mi,m 2 J = , t~- 



\\K\\ 



{2.1 



where 



K 2 + m 2 -m 2 2 _ v /A[K,m 1 ,m 2 \ 
a[K } m l ;m 2 ] = —. , p[K; m 1} m 2 \ = . (2.9) 



K 2 



K 2 



Notice that when mi = m 2 = we have a — (3 — 1, thus reproducing the massless case. Finally 
the original £ can be parameterized as 



tP x ~ x + zK 



K 2 



\\K\X 



XX 



K 2 



K )+a- 



XX 



K 2 



K 



(2.10) 



3 If E < 0, we need 24 
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2.2 Spinor integration method for the physical phase space integration 

Now, we explore how to apply the spinor integration method for massive cuts to the phase space 
integration. Just like the massless case, when there are only two outgoing particles, spinor inte- 
gration method can be applied directly without any modification. To see explicitly, just write the 
phase space of the cross section: 



'»= II /^-»!)M'*'(i'-Ey 

f=l,2 J K ' f=l,2 

~ J d*L x 5 + (L\ - m\)5 + ({K - Lx) 2 - m 2 ), (2.11) 

which is exactly the same (namely the show-up of two 5-functions) as the spinor integration method 
given in Eq. (|2.1| ). 

Thing will be different when n = 3, where the physical phase space is given by: 

J 3 = II / 0f3 6+ ^ - m ^)^ K - E ^ L 3) 

= / ^f 5+{Ll ~ ml) J ^ + (Ll-rnl)6^(K - L 3 - L 2 ) 2 -ml)f(L 2 ,L 3 ) 

= j ^ + {Ll-ml)f{Lz). (2.12) 

The problem we meet here is just the same as the massless case. The integration over L 2 with two 8- 
functions can be performed by the spinor integration method. However, there is only one ^-function 
in the integration over L 3 . In order to apply the spinor method recursively and continuously, we 
need insert one more 5-function like the Faddeev-Popov method. 

Similarly to the massless case [|15|, we consider the following integration 



I x = J dx5{{xK - L 3 f -ml) = J dx5(x 2 K 2 - x(2K ■ L 3 ) + L\ - 



where the 5-function has two solutions 



(2K ■ L 3 ) ± ^{2K-L 3 f-AK\Ll- 



2\ 



Xi = v- - ^ -j* (2 13) 



Using the on-shell condition L\ = m 3 , it reduces to 



2K-L 3 , 
x- = 0, x + = — ^ . (2.14) 
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We find that x_ = is always a root. However, when x = 0, we have 8{L\ — m§) which does not 
give an independent ^-function. So x_ = should be excluded from our consideration. For another 
root x + , from (K — L 3 ) 2 = K 2 + m| — 2K ■ L 3 > (m x + m 2 ) 2 , we have 

2K ■ L 3 < K 2 + m 2 - (mi + m 2 ) 2 , (2.15) 

which gives the upper bound of x + . For the lower bound, considering the center-of-mass frame 
where K = (E, 0,0,0), L 3 = (E 3 ,p,0,0) with E\ - p 2 = mj, we have 2K ■ L 3 > 2Em 3 , i.e., 
(2K ■ L 3 ) 2 > AK 2 m 2 . 

Putting all consideration together we have 



rxi rx 1 

I x = / dxS((xK — L 3 ) 2 — m 2 ) = / dx- 

J In J In 



2x+K 2 - (2K ■ L 3 



I x J x 

1 

W{2K.L 3 ) 2 -AK 2 {L 2 -mf ) 



(2.16) 



where 



4m| K 2 + m 2 — m 2 



x [K,m 3 ] = Xi[K,m 3 ;m t ] = ^| = a/ x q [K, m 3 } 2 + A[iT; m 3 , m t } 2 ,(2.17) 



where A[iT; m 3 , m t ] = a/A [fT, m 3 , m t }/K 2 with m ( = mi + m 2 . Using i^ 2 > (mi + m 2 + m, 3 ) 2 , it is 
easy to see that A[K; m 3 , m t } 2 > and thus x± > x . 
Now Eq. ( [2.1 2j ) can be written as 

/ d 4 L 3 6+(Ll - m 2 3 )\^/(2K ■ L 3 ) 2 - AK 2 {L 2 - m 2 )| H dx6((xK - L 3 ) 2 - m 2 )/(L 3 ). 

J J xn 



(2tt) 3 

Decomposing L 3 = £ + zK with £ 2 = 0, then 

I 3 = — ^- / cfed 4 £5 + (£ 2 )(2£ • K)5 + (z 2 K 2 + 2zK ■ £ - m 2 )(2L 3 • K) j * dx5 + (x 2 K 2 - 2xK ■ L 3 )f(£) 

y^ 71 ) J Jx 

= — ^ / rfxxK 2 f dz(x - 2z)K 2 5 + (z(x - z)K 2 - m 2 ) 

x J d 4 £5 + {£ 2 )5 + {x{x-2z)K 2 -2xK -£)J{£). (2.18) 
One by-product of the above procedure is 

2K-L 3 = xK 2 . (2.19) 
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By solving the 5-function 5 + (z(x — z)K 2 — m 2 ) and the similar kinematical discussion as in Section 
2.1, we get 



xK 2 - ^K 2 (x 2 K 2 -Am 2 



2if 2 



x — 2z = \l x 2 — x\ . 



(2.20) 



Continue the evaluation as 



(2vr) 3 



(2tt) 3 



(2tt) 3 



j-i 



X 
XI 



dxxK 2 / {X\d\) X\dX tdU + (x(x - 2z)K 2 - xt (X\K\X )f{X,X,t) 



XQ 



dxxK 2 I (X\dX) 



XldX 



J-l 



dx{K 2 f\ x 2 -x 2 I (X\dX) 



x(X\K\X 



X\dX 



X\K\X 



K 2 j. 



i~t " 2 / >i 2 



AIATIA 



1 2 ' 



AIKIA 



(2.21) 



where c = 7r/2 is related to the Jacobi of changing integration variables and the way we have taken 
the residues. 

Finally we arrive 

d A U 



*5 + (Ll-mj)f(L 3 ) 



(2tt)< 



7T 



2(2tt) 



dx(K 2 ) 2 J x 2 - x 2 / (X\dX) X\dX 



/(A, A, t) 



X\K\X 



2 ■ 



I = K y X2 J*0 P^22) 



X\K\X 



This is our key setup in this paper. Notice that when mi = m 2 = m 3 = 0, Xq = 0, so Eq. ( [2.22 ) 
reduces to 



d 4 L 3 



7T 



2(2vr) 3 



S + (Lj)f(L 3 



dx(K 2 ) 2 x / (A|dA) X\dX 



/(A, X,t) 



K 2 x 



X\K\X 



2 ■ 



X\K\X 



(2.23) 



which is the familiar massless case presented in [15|. 

In the end, let us give a remark. The integration region of x G [xq, X\] depends on the dynamical 
momentum K as well as mass parameters m 3 and m tota i. Because this, the roles of m% and m to tai 
are not obviously symmetric. If we define 



x = \ Xq + A 2 u 2 



(2.24) 
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the integration region of u will be u G [0, 1] which does not depend on external momenta and masses 
anymore. Under this transformation we have 

A 3 u 2 



Ti 



2(2tt) 



du(K 



2\2 



'x\ + A 2 u 2 



(X\d\) \\d\ 



/(A, A,*) 



K 2 Au 



\\K\X 



2 ■ 



\\K\X 



(2.25) 



and 



L3 = J^ AA 



y^+AV Au 



X\K\\ 



K. 



(2.26) 



This transformation will become even simpler when = where we get a just linear transformation 
x = Au. Although Eq. ( |2.25| ) may look simpler, however by some calculations we find that Eq. 



Q2.25Q is, in general, not better than Eq. ( |2.22 ) and readers can use anyone they like. In the later 
part of this paper, we will use the form of Eq. ( 2.22j ) . 



3. Simple examples 

In this section, we present some very simple examples to demonstrate our method, especially the 
integration region of x. We denote the physical phase space integration of n outgoing particles as 
In° rm {f'i where s stands for the spinor method and m the momentum method. The K is the 
sum of momenta of these n particles and / is a general function. 

3.1 The pure phase space integration with two outgoing particles 

This integral can be performed directly by the spinor method as we have analyzed in last section. 
Spinor integration method: The integration is given by 



d A L 2 d 4 L 



1 5 + {L\ - m z 2 )S + (L( - m{)(2Txf5\K -L 2 -U 



(2tt) 3 (2tt) 3 

= (2^5 / d ' L ^^\ ~ <)t + ((K - L,f - ml). 
According to Eq.( [2.7|) , one gets 

1 - 2z)K 2 + m\ - ml 1 vr 



(3.1) 
(3.2) 



2(2tt) 



(A A) 



A A 



\\K\X 



(2vr) 2 2 







1 y/(K 2 -ml- m 2 2 ) 2 - Am\ 



2 mo 



2(2tt) 1 



K 2 



(3.3) 



S 



which is obviously symmetric between m\ 1 m 2 . 

Momentum integration method: It is given by 



dL\ 



(2vr)V(K-L 2 -L 1 ). 



(2tt) 3 2£i (2tx) 3 2E 2 

Taking the center-of-mass frame, where K = (E, 0, 0, 0) and L\ = (Ei, ki, 0, 0), yields 



I™(1;K) 



(27T) 



k\dk\ 
2Eo 



dVt5 + (E 2 - 2EE X +m\- m\) 



1 7T 

(2tt) 2 2 



0- 



(3.4) 



(3.5) 



3.2 The pure phase space integration with three outgoing particles 



Spinor integration method: From Eq. fl2.12p with the result I 2 (l; K) in previous subsection we 
have 



d A L, 
(2n) : 



■5 + (Ll-ml)I°(l;K-L 3 
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2(2tt) 



j-i 



dx(K 2 fJx 2 -xl / (A|dA) \\d\ 



•t'li 



Ii(l;K-L 3 ) 



\\K\X 



2 ' 



(3.6) 
(3.7) 



where If (1; K — L 3 ) depends on (K — L 3 ) 2 only. But using L\ — ml = and (xK — L 3 ) 2 — ml = 0, 
we can find 



(K - L 3 f 

which does not depend on A, A at all. Thus 



[1 - x)K 2 +m 2 3 , 



7i 



2(2tt) ; 



xi 



.'■() 



dxK 2 \/x 2 -x 2 I°(1;K-L 



(3.8) 



(3.9) 



where xx\ is given by ( J2.17 ) . This expression is obviously symmetric between mi,m 2 , but not for 
m 3 . However, it is easy to check numerically that the final result is indeed symmetric among all 
the masse parameters. 

Momentum integration method: The integration is 

d 3 L x d 3 L 2 d 3 L 3 



(2tt) 3 2E 1 (2tt) 3 2E 2 {2tt) 3 2E 3 
d 3 L 3 



(27r) 4 5 4 (K-L 1 -L 2 -L 3 ) 



(27T) 3 2£; 



■I?(1;K-L S ). 



(3.10) 
(3.11) 



In the center-of-mass frame, K = (E, 0, 0, 0), L 3 = (E 3 ,p, 0, 0) with E 2 -p 2 = ml, thus (E-E 3 ) 2 - 
p 2 > (mi + m 2) 2 , i-e., E 3 < (E 2 + ml — (mi + m 2 ) 2 /2E). Namely, the integration region of E 3 is 

E 2 + m 2 3 — (mi + m 2 ) 2 



m 3 < E 3 < 



2E 



(3.12) 
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Using this we have 

4 m (l; K ) = J^yf dE^El-miI™{\- K — L 3 ). (3.13) 
In order to show this is identical to Eq. ( |3.9| ), we can make a transformation 2E 3 /E — > x. Then 



£ 2 + m§ - (mi + m 2 ) 2 . . 

™s < £ 3 < > xa<x <xi, 



2EyJEl - m 2 if 2 -^ 2 - xg. (3.14) 

It's obvious to see that 7^(1; if) = J|(l, if). 

3.3 The pure phase space integration with four outgoing particles 

Here we will only present the expression using spinor method. The pure phase space is 

m\ K ) = j f^l 5+ ( L 4 - om, k-l a ). (3.i5) 

Using the recursive method we get 

„( 4 ) 



.0) 
-1 



dx^\K - L 4 ) 2 V /(x( 3 )) 2 - {xt ] f I s 2 {\- K-Ls- U). (3.16) 





where naively we have following boundary values 



(3) _ / 4m| (3) _ {K- Uf + ml - (m, + m 2 ) 2 



However, similarly to Eq.( |3.8| ), we can find that 

{K-Uf = (l-x^)K 2 + ml 
(K - U - L 3 ) 2 = (1 - x (3) )(if - L 4 ) 2 + m 2 = (1 - x (3) )(l - x (4) )if 2 + (1 - x^)ml + m 2 (3.18) 

Thus we have 



(4) / 4m 2 {4) _ K 2 + ml-( mi +m 2 + m 3 ) 2 



(3) / 4m 2 (3) (1 - x (4) )if 2 + m 2 + m 2 - (mi + m 2 ) 2 

X ° M/ (l-x( 4 ))K 2 + m 2 ' Xl (l-x( 4 ))K 2 + m 2 1 j 
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Putting ( 3.18|) and ( |3.19|) into ( |3.16 ), we get the analytic expression for the pure phase space of 
four arbitrary massive particles 



,.(4) 



(4) 




dx {4 \ K 2 {(x^) 2 K 2 -Am 2 



„(3) 



dx^ 



(3) 




64(2tt) 5 

[((1 - x^)K 2 + m 2 )((x^) 2 ((l - x^)K 2 + m 2 ) - 4m§)]a 



x 



(1 - x( 3 ))(l - x^)K 2 + (1 - x^)m\ + m\ 



x[((l -x (3) )(l -x (4) )K 2 + (1 - x {z) )m\ + m\ - m\ - ml) 2 - Am\m 2 ^ . (3.20) 

The expression is not obviously symmetric among (mi,m2,m 3 ,m4 : ) by our choice of the order of 
integrations. However, it is easy to check by the numerical method that the final result is indeed 
symmetric among all masses. 

3.4 The phase space integration of 3-outgoing particles with / = (2L X • L 2 )(2L 1 ■ L 3 ) 

Here we calculate a relatively complicate example with / = (2L% ■ L 2 )(2Li ■ L 3 ). 
Spinor integration method: The integration can be directly written as 

d 4 U 



J3 5 + (Ll-m 2 3 )m;K'), K' = K — L 3 



(2tt)' 



where using the momentum conservation, / can be written as 

/ = (2Lx ■ K' - 2m\)(2L 1 ■ (K - K')). 
Using the Eq. (|2.10|) with K replaced by K', we can simplify / further as 

'\\K\X 



(3.21) 



(3.22) 



/ = (aK' 2 - 2ml) P K 



/2 



where 



a 



K' 2 + m{ — m 2 



\\K'\\ 



+ (a - f3)K' -K -aK 



^ A[K',m 1 ,m 2 \ 
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(3.23) 



K' 2 ' ' K' 2 
Now we calculate K') using the simplified version /. It is given by 



W;K') 
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a(3{K' 2 -mf- m 2 ){K' ■ K — K' 



2(2vr) 2 

When we put it back into If (/; K), we need to know that 

K' 2 

K' -K = K 2 - zK 2 



(3.24) 



(3.25) 



1 - x)K 2 + ml 



1 



t(X\K\X 



K 2 - -xK 2 , 
2 
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where we have used the relation t = ^h^jjp This means K') does't contain explicitly A and 
A for the spinor integration over L 3 . So we get immediately 

W; K) = £ dxK 2 ^/x 2 -x 2 aP((l - x)K 2 + m \- m \- m^-xK 2 - m\). (3.26) 



where xq,Xi and a, (3 are given by ( |2.17| ) and ( |3.24j) respectively. 



Momentum integration method: The integration is 

d 3 L 3 



I?U\K) = I J^4e/2V;K'). (3.27) 



I™(f] K') can be calculated as follows: 

Choose a center-of-mass frame, such that K' = (E', 0,0,0), L\ = (Ei, 0, 0, k\) with E\ — k\ = m\ 
and K = (E,p). The angle between p and Li is 9. Then 

x(2E 1 E' - 2m\)(2(EE l - ypk x ) - 2E 1 E') 
= — I hdEiS + (E' 2 - 2E'E 1 + m\- m\){2E x E' - 2m\){2EE 1 - 2£i£') 

27T y 

_ 1 E' 2 + mf — m\ ^A[E>,m umj ] {E , 2 _ ^ _ ^ _ ^ 



In E' E' 2 



This is identical to if') in the center-of-mass frame. To calculate I™(f;K) simply, we need 
to choose the center-of-mass frame of K which is not the one we have used for I™{f\ K'), thus we 
need to write i™(/; K') as the Lorentz- invariant form, which is not so straightforward sometimes. 

Here we use the Lorentz invariant form I^fjK') given by the spinor method to go further. 
Taking the center-of-mass frame where K = (E, 0, 0, 0), then 



m, , x 1 f ,^ l~, E* - 2EE 3 + mi + mi - mi r— T 

x(E 2 - 2EE 3 + m\-m\- m 2 2 )(EE 3 - m\). (3.30) 

It is equal to i|(/; K), which can be easily checked by making a transform 2E 3 /E — y x as in Section 
3.2. 
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4. Practical applications 



In our previous section we have done some simple examples. However, these examples do not involve 
the real amplitudes. In this section we will discuss the phase space integration of two simple real 
physics progresses with three out-going particles. These two examples are presented in the following 
two references: [|B| and [P0fl . 

4.1 Z° decays into lepton pairs and spin-0 bosons 

This example discuss the decay reaction 

Z° ->■ 1+rH, (4.1) 

where / stands for the electron or muon with mi = and H for the Higgs with mn ^ 0. The invariant 
matrix element squared is given by Eq.(2.10a) in [IS]. According to the Glashow-Weinberg-Salam 
model (Eq.(3.1) and Eq.(3.2) in the matrix element squared can be written as 



pol 



rM 2 9 



where k is the total momentum of Z°, Mz the mass of Z° and 



9 



- sin 2 6 



w j 



9 



Q = l+ + L 



cos W A 4 cos 6 W M z cos 9w 

We evaluate this by first evaluating the phase space integration over Z + and Ih- From Eq. 
can easily get 

(k-i-y 



pol 



71 



2(2tt) 



(A dX) 



X dX 



n 2 



X ■ 



(t ( A|Z_|A 

(A dX) 



C 7T 



2(2tt 



A dX 



X\k-l-\X 
(t{X\l_\X M 2 z + 2{k-lJ)t{X\k\X ) 

{{k-l-Y-Mlf 



AlPilA 



A|P 2 |A 



X\R\X 



where 



Pi 

R = l_ + 



(k-l-)- 
2(k-l 



(k-lA 2 -Ml x 



Po = k-L 



Ml 



-k. 



(4.2) 



we 



(4.3) 
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Introducing one Feynman parameter we can continue to 



Ja(l>-1- 

pol 



C IT 



2(2tt)2 J Q 



dy I (XdX) XdX {{k-l_f -M 2 H f 



C 71 



where 



2P-R = 2k-L 



2(2vr)2 J q 
2(k-l 



1 dy{{k-l-f-Ml^ y2P ' R 



2y(X\R\X 



X\P\X 



l 3 



(P 



2\2 



Ml 



(sy + u + k ) 



P 2 = sy + u, s = —r^k-l-f - {k 2 -M 2 H )(2k ■/_)), u = k 2 - 2k ■ 1. 



Now we put Q4~4" ) into 



js(I>) 



d 4 L 

(2tt) 



;5 + (i 2 _mj2;k-i-). 



pol 



To continue, we exchange the integration order of J dy and J dH-. Using 2k ■ £_ = t (X\k\X 
and performing the spinor integration, we finally arrive 



pol 



16(2tt) 







2 \2 



dx(k 2 x) 2 / - rx - Mfj) 



y(l + ^(sy + u + k 2 )) 



(sy + u) 



where 



k 2 x 
M 2 



{k 2 {l-x)-M 2 H ), u = k 2 {l-x). 



We can perform the integral over y to yield 

k 2 -Mfj 

hC£;k) 



pol 



16(2tt)3 J 

(1 - x)M 2 {k 2 + Ml 



dxMl {k 2 x + M 2 Z ) In 1 



x(k 2 (l — x) — Mjj) 
(l-^)Ml 



k 2 x(k 2 (x- 1) + Af2) 
+ V . ,o — -k 2 -Ml 



Ml 



(1 - x)(k 2 x - Ml) -xM 2 H 

which can be integrated further to get analytic expression if one wants. Notice that k 2 = M z 
fl4.7|) can be simplified further as 

„2 . / A^tr 



Js(5>) 



2M I / 2 ,2> 

z (a 2 + b 2 



g 



pol 



384(vr) 



Ml cos 2 W y 



i_ M i 



dx I (x+ 1) In ■ 



l + ar + x\ 



Ml 2 ) 



1 — X 
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2(1 -x) 



+ X + x( 



l + x 2 + x(j^-2 



H 

M 2 



-1: 



(4.8) 



Multiplied by the normalization factor (2Mz)~ l we have omitted in the calculation of the cross 
section, the integrand of Eq. (f!8|) is just Eq. (3.4) in [|19|] by verifying x = X- and 5 = Mh/Mz in 
the center-of-mass frame. 

4.2 The production of Higgs bosons in pp collisions 

For the second real example, we consider the quark-antiquark-annihilation mechanism qq — > ffH° 
in |p0[ . The corresponding cross section and the matrix element squared are respectively given by 
Eq.(2.2) and Eq.(2.3) in PD| . We write the cross section as 



I°(H^q» u ; Q) = c J d 4 k5 + (k 2 - m)) J d 4 k5 + (k 2 - m)) J d 4 h5 + (h 2 - m 2 H )5 4 (Q -k-k- h)H^q^ 



c J d 4 k5 + {k 2 - m 2 f )P 2 (H^q^- Q'), Q' = Q-k, 



(4.9) 



where k(k) and h are respectively the momentum of the heavy quark(antiquark) /(/) and the 
Higgs, Q, q, q the total momentum and momenta of particles q, q. We have absorbed all the common 
constant factors including 7r-factor into c. H^q^ is given by 



E^q 



fJ,U 



32 



(2h ■ k + m 2 H )(2h ■ k + m 2 H ) 



Q 2 (Q-hf 



(Am 2 - m 2 H )Q 2 



1 + 



(2h ■ k + m 2 H )(2h ■ k + m 2 H ) 



(Q 2 + m 2 H - 4m 



2Q ■ h(4m 2 - m 2 H/ 
(2h-k + m 2 H ) 



Q 2 

T'"< 



m 2 f — 2k ■ q k ■ q 



+ 



2Q ■ h(4m 2 - m 2 H ) 



(Q 2 + m 2 H - Ami) + 



Notice that 



fJ (2h ■ k + m 2 H 
- (Q 2 + m\ - Am 2 f )[2k ■ q k ■ q + 2k ■ q k ■ q - Q 2 k ■ k] 



2h-k + m 2 H = Q 2 -2Q-k 
2h ■ k + m 2 H = Q 2 - 2Q ■ k 

2Q ■ h = (Q 2 - 2Q ■ k) + (Q 2 -2Q-k) 



Q 2 

— m 2 — 2k ■ q k ■ q 



(4.10) 



2k ■ k = (Q 2 



m 



H 



2mi 



(Q 2 -2Q-k)-(Q 2 -2Q-k). 
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To simplify the calculation, we rearrange H^q^ as 
1 „o \ Q 2 - 2Q ■ k Q 2 - 2Q -k 



H nu g »u = 32 ^ _qi 



Q 2 -2Q-k 



2\2 



+ -(Amj-mj I )(Q 2 ) 



Q- 



Q 2 -2Q ■ k 
1 

+ 



(Q 2 -2Q-k) 2 (Q 2 -2Q-k)(Q 2 -2Q-k) (Q 2 - 2Q ■ kf 



+ 



Am 2 — m 2 H 



l(Q 2 - 2Q ■ k)(Q 2 -2Q-k) (Q 2 - 2Q ■ k) 



[Q 2 m) -2k-q2k-q] 



Q 2 + m 2 H — Am 2 , 



(Q 2 - 2Q ■ k)(Q 2 - 2Q ■ k) 
(Q 2 + m 2 H - Am 2 )Q 2 



\2k-q2k-q- -Q 2 {{Q 2 + m 2 H - 2m 2 ) - {Q 2 - 2Q ■ k))} 
I 



2(Q 2 -2Q ■ k) 



(4.11) 



where we have used the symmetry between k and k. 

Now we can start the calculation. First we will perform the phase space integration over k and 
h. Then we perform the left k integration. 



4.2.1 The integration I°{H^q^ u -Q') 

First we simplify the input according to Eq.( |2.1U| ), i.e. 

Q' 2 



k 



Q 2 -2Q-k 



X\Q'\X 
Q' 2 



' \ XX T^—y ) + Q,2 Q 



Q r 



X\Q'\X 



AlPilA 



(4.12) 



with 



Pi 



Q' 



By checking Eq.( 4.11| ), we find that there are four nontrivial integrations should be attacked. We 
do them one by one. 

Type I: f 1 = l/(Q 2 -2Q-k) 

The integration is 



(A dX) 



X dX 



(3ai 



AKQ'|A 



AIPJA 
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dy / (A dX) 



X dX 



/3cti 



X\yP 1 + (l-y)a 1 Q'\X 



1 , 4/^Q' 2 
dy- 



where 



o y -y 2 f3 2 ^ + Aa 2 (Q> 2 ) 2 ' 



Z = (2Q-Q') 2 -4Q 2 Q' 2 . 



Type II: f 2 = Q 2 -2Q-k 
The integral is 



IltftQf) = / (AdA) 



A rfA 



PQ' 2 Q' 2 



X\Q'\X 



X\Q'\X 



A|Pi|A 



= /3Q' • Px = (5 ai Q' 2 . 
Type III:/ 3 = l/(g 2 -2Q-fc) 2 



Wastf) = f (A dA> 

= y (a dx) 



X dX 



X dX 



12 



1 



A|Q'|A 

NQ' 2 

A|Pi|A 



2 (Q 2 - 2Q ■ kf 



4/3 



-/3 2 £ + 4a 2 (Q 



Type IV: / 4 = 2k ■ q/{Q 2 - 2Q ■ k) 
Using 



2k ■ q -- 
we get 

Ws;Q') = [ (AdA> 



/2 



A|Q'|A 



A rfA 



(A|P 2 |A 


, p 2 = / 


PQ' 2 


(A|P 2 |A 


(a|Q'|A 


2 (a|Pi|a 
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2PalQ'\l-y)(\\P 2 \\ 



l 3 



(4.20) 



dy J (A dX) X dX 

X\yP 1 + {l-y)a 1 Q'\X 

1 ^IQ'\1 - l/)[(-l//3 2 (l - jgK aia)2g ■ g - + a 1 a)2fc ■ g] 

* (-Z/ 2 /3 2 4j^ + «?Q' 2 ) 2 ' 

Now substituting above four types of integrations into I^H^q^; Q') and with some algebraic 
manipulation, we can get 

1 / 2 «(iFV"; Of) = [W - 20 ■ t) + I(^M-^L _ if - W + m\ 



32 



L 4 



Q 2 -2Q ■ k Q 2 -2Q-k 



-Ami) + -Q 



1 2 (Q 2 + - 4m 2 ) (Q 2 + m 2 H -2m 



Q 2 -2Q-k 



dy- 



4/3a 1 Q /: 



-l/ 2 /3 2 S + 4a 2 (g 



(g 2 + m 2 f -4m 2 f )(2fc-g) Z" 1 



^/3a 2 g ,2 (l-y) 



g 2 - 2g • k 

[(-Z//3 2 (l - + «i«)(2g • g) - (y/3 2 ^# + a x a)(2* • g)] 



-y^^ + ^g' 2 ) 2 



(3 (Am 2 — m 2 H ) 
-/3 2 S + 4a 2 (g /2 ) 2 



[(g 2 ) 2 + 4(g 2 m 2 -2A;-g2A;-g)] + 



(4m 2 - m 2 ,)/3(g 
4(g 2 - 2g • fc): 



2\2 



g 2 /3 ai g /2 - 2g 2 /3(g 2 + m 2 H - Am 2 
+ 4(g 2 -2Q ■ k) 

4.2.2 The integration I£(f;Q') 

Now we do the left k- integration using fl2.22| ) with 



+ 



(4.22) 



and following relations 



x = 

Q-k 
2k ■ q 

2k ■ q 



'Am 2 Q 2 + m 2 — (mj + mnf 

g r ' Xl= g 2 



(4.23) 



x 



Q\ (Q - kf = (1 - x)Q 2 + m 2 



(x-22i)(A|g|A 



^A|g|A 

(x-2z 1 )(A|g|A 



A|g|A 



+ 2z x q ■ q 



2ziQ ■ q 
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From I 2 (H^ 1 'q^ 1 Q'), the terms containing Q • k and (Q — A;) 2 do not depend on k, thus can be 
done easily just as in the example of the pure phase space integration with three outgoing particles. 
Then we need only perform the following two types of nontrivial integrations: 



mq-k;Q) 



dx(Q 2 ) 2 \ x 2 -x 2 / (A|dA) X\dX 



1 (x-2z 1 )(X\q\X 

Zt~2 ( 



X\Q\X 



X\Q\X 



Q 2 + 2z x Q ■ q) 



Xl 



xo 

XI 



dx(Q 2 ) 2 Jx 2 -x 2 [(x~2 Zl )^ + 



Q ■ q 2ziQ ■ q 



Q 2 



Q 2 



/ dxQ 2 \J x 2 — x\xQ ■ q, 



(4.24) 



I°(2k- q2k- q;Q) 



dx(Q 2 ) 2 J x 2 - x 2 / (X\dX) X\dX 



XI 



{x-2z 1 ) 2 (x\q\X 




(X\Q\X 


2 



+2 Zl Q ■ q- 



{x-2z 1 )(X\q\X 



X\Q\X 

dx{Q 2 ) 2 \/x 2 -x 2 ((x-2z 1 ) 2 
Q-q 



(Q 2 Y + 2z 1 Q-q- 



Q 2 + 2 Zl Q ■ q2 Zl Q ■ q) 



X\Q\X 

{x-2z x ){X\q\X 



X\Q\X 



Q 2 



{2Q-q){2Q-q)-Q 2 {q-q) 

3Q 2 

Q ■ q 2ziQ ■ q2ziQ ■ q, 



+2z x Q ■ q(x - Zzi)^ + 2^iQ ■ q(x - 2z{)^- 4 



Q 2 



X'() 



dxQ 2 J x 2 — x"l 



(2Q • q)(2Q • qJ {x2Q2 _ _ q_q {x2Q2 _ . ) 



3Q 2 " " J/ 3 

4.2.3 The final result 

Substituting I S 2 {H» V q» v ; Q'), l£(2k -q2k-q;Q) and J 3 s (2g ■ k; Q) into Eq.QOQ yields 



dx \ x 2 — Xr 



X'O 



dy 



P{2 -2a + ax) 



°{2(l-x) J -y 2 (3 2 (x 2 - x 2 ) + (2 - 2a + ax) 



(Q 2 + m 2 H - 4m 2 ,) f 1 



dy- 



Af3(l-y){2-2a + ax) 



Q 2 (l-x) J y (-y 2 f3 2 (x 2 -x 2 ) + (2-2a + ax) 2 ) 2 
f3(4m 2 - m 2 H )/Q 2 



To 



-@ 2 {x 2 - x 2 ) + {2-2a + ax) 2 



T 3 + T< 



4 r, 



(4.26) 
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7 3 7(32c) ([GeV] 4 ) 




Figure 1: 3^/3 as a function of -y/s for the pp collision. The dashed and continuous curves respectively 
represent two sets of parameters: m# = 10 GeV, mf = 4.5 GeV and ran = 30 GeV, nif = 35 GeV. 



where 



T x = (Q 2 x + m 2 H - Am)) 2 + (Q 2 + m 2 H - Am)) 2 + A(Q 2 + m 2 H - Amj)mj + 2Q 2 (6mj - m 2 H ) 

+ 3 (q ■ q(x 2 Q 2 ~ 4mJ) - (2Q • q)(2Q ■ q)(x 2 - -^)), 
2m 2 

T 2 = (yP 2 (x - —f) + a{2 -2a + ax))x{Q ■ q)(Q ■ q) 
Q 



+ -(yf3\2 -x) + a(2 -2a + ax))(q ■ q(x 2 Q 2 - Am 2 ) - (2Q ■ q)(2Q ■ q)(x 2 - -^)), 
T 3 = (Q 2 ) 2 + AQ 2 m 2 + % ■ q(x 2 Q 2 - Am 2 ) - (2Q • q)(2Q ■ q)(x 2 - )), 



_ (3(Am 2 f - m 2 H )Q 2 (3(2 - 2a + ax)(Q 2 ) 2 - A(Q 2 + m\ - Am 2 )Q 2 (3 1 
Ti ~ 4(1 - xf + 8(1^) + Y {Q ] ■ (427) 

The corresponding parameters are 

= (1 - x)Q 2 + 2m 2 - m 2 H = [((1 - x)Q 2 - m 2 H ) 2 - Am)m 2 H \\ 
a ~ (l-x)Q 2 + m) ' 1 (l-x)Q 2 + m) ' { ] 

In Fig. [I], we display g^r-Tf versus the cm. energy y/s of the pp by the numerical method. 
Notice that the displayed 3^/3 is not the real cross section since the dynamical factor c given in the 
original reference has not been included and the real cross section is 7| (so the decay behavior of 
1% at high energy can not be observed from this Figure). Here, we emphasize two points. First, by 
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our spinor method, almost all calculations have been reduced to reading out the residues of poles 
and making some algebraic manipulations. Thus although the analytic expression looks long, the 
calculation is kindly trivial. 

Second, we can take appropriate integration order to simplify the process according to the 
structure of the integrand. Usually we should first perform the integrations over those variables, 
with respect to which the structure of the integrand is relatively simple. In this example, we have 
leave k as the last integration variable 4 . This is because the integrand Eq. (|4.11|) does not contain 
h explicitly. 

Notice that different integration ordering, i.e., integrating over p\ first and then p 2 or integrating 
over P2 first and then pi, will in general give different-looking expressions. For example, in the 
expression of Q3.20 ), we have fixed arbitrarily the ordering m 1 ,m 2 ,m3,m 4 . Different ordering will 
end up with different integration regions although the final result should be the same. Furthermore, 
if we have left one particle un-integrated while others have been integrated, then we will get the 
corresponding differential cross section for this particle. Thus different integration ordering will give 
different differential cross sections for different particles. 



5. Conclusion 

Originating from the application of the spinor method to the massless case, in this paper, we have 
established the framework to process the massive case. From the examples presented in the paper, 
the advantages of our method is further manifested. 

First, the manifestly Lorentz invariant form of the result in each step is gotten naturally. 
This ensures that the recursive method can be applied conveniently especially when the number 
of outgoing particles is large. In this process, we don't need take any specified reference frame as 
when using the momentum integration method. 

Second, the integration regions can be written straightforward according to Eq.(2.17), while 



with the momentum integration method, one has to pursue exhaustively to specify those of many 
variables (for example, angles and module variables). Note that in our method, though for the 
massive case the region is not so simple as the massless case, it is only the functions of mass and 
energy. 

Finally, the salient point is that the constrained three-dimensional momentum space integration 
is reduced to an one-dimensional integration, plus possible Feynman integrations. However, in this 
so large simplification, we just pay a little extra price, namely the integration over A and A which 
can be obtained by reading out residues of corresponding poles. 



4 Because the symmetry of k and k, it's the same if we leave k as the last integration variable 
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In this paper, our new method have shown out the value of practical calculations. As we have 
mentioned in the introduction, our method provides compact analytic expressions for cross section. 
Thus we can investigate the analytic structure using these expressions. We think it is an interesting 
direction. Also, in this paper we have just touched the tree-level result. It is our goal to combine 
these analytic expressions with one-loop result to see if we can improve current numerical NLO 
algorithm, especially the IR singularity substraction. A regularization scheme is mandatory and we 
need consider the general D-dimensional case. All these questions will be our future projects. 
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